Study of light-cone distribution amplitudes for p-wave 

heavy mesons 

Chien-Wen Hwang* 

Department of Physics, National Kaohsiung Normal University, 
Kaohsiung, Taiwan 824, Republic of China 

Abstract 

In this paper, a study of light-cone distribution amplitudes for p-wave heavy mesons is presented 
in both general and heavy quark frameworks. Within the light-front approach, the leading twist 
light-cone distribution amplitudes, <Pm(u) and their relevant decay constants of heavy scalar, axial- 
vector and tensor mesons, /m, are formulated. The relations of some decay constants can be 
simplified when the heavy quark limit is taken into account. After fixing the parameters which 



appear in a Gaussian wave function, the corresponding decay constants are calculated and compared 
with those of other theoretical approaches. The curves and the first six £ moments of 4>m(u) are 
plotted and estimated. These results all endorse the requirements of heavy quark symmetry. 
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I. INTRODUCTION 



Light-cone distribution amplitudes (LCDAs) of hadrons are key ingredients in the descrip- 
tion of various exclusive processes of quantum chromodynamics (QCD), and their role can be 
analogous to those of parton distributions in inclusive processes. In terms of Bethe-Salpeter 
wave functions (p(ui,ki±), LCDAs are defined by retaining the momentum fractions 

Ui and integrating out the transverse momenta ki± They provide essential information 
on the nonperturbative structure of the hadron for QCD treatment of exclusive reactions. 
Specifically, the leading twist LCDAs describe the probability amplitudes to find the hadron 
in a Fock state with the minimum number of constituents. In the literature, there have been 

m, 

Nabmbu- 



many nonperturba' ;ive_a.] >] m >aehes to cstnua.tr LCDAs. such as the QCD sum rules 



lattice calculation 7, 8 
Jona-Lasinio model 111 



dealt with LCDAs of pseudoscalar 3l. l8Ml4i1 . vector 

n 

tensor [6( 



chiral quark model from the instanton vacuum 
and the li ght- front quark model 



m 



13 



13 



151 ] . These studies have 
14j , axial- vector jl], [f], 2] , and 



mesons. 

The fact that .B-physics exclusive processes are under investigation in BABAR, Belle, and 
LHC experiments also urges the detailed study of hadronic LCDAs. Recently, many p-wave 



heavy mesons were observed and confirmed. They include Dq°, D\, D^, Dfx, D* 2 16M20| 
and B®, B^ , B° sl) B*® |2l|-|23|. The present paper is devoted to the study of leading twist 
LCDAs of p-wave heavy mesons which include the scalar, axial-vector, and tensor mesons. 
We hope a thorough understanding of their properties, such as LCDAs which are universal 
nonperturbative objects, will be of great benefit when analyzing the hard exclusive processes 
with heavy meson production. 

In the past decade, the most significant progress made in the QCD description of hadronic 
physics was, perhaps, in the avenue of heavy quark dynamics. The analysis of heavy hadron 
structures has been tremendously simplified by the heavy quark symmetry (HQS) proposed 
by Isgur and Wise 24], |25[ , and the heavy quark effective theory (HQET) developed from 
QCD in terms of 1/rriQ expansion HQET has provided a systematic framework for 

studying symmetry breaking 1/mg corrections (for a review, see Ref. [291). Moreover, in 
terms of heavy quark expansion, HQET offered a new framework for the systematic study 



of the inclusive decays of heavy mesons 



3CH33|. However, the general properties of heavy 



hadrons, namely, their decay constants, transition form factors, structure functions, etc, are 
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still incalculable within QCD, even in the infinite quark-mass limit with the utilization of 
HQS and HQET. Hence, although HQS and HQET have simplified heavy quark dynamics, 
a complete first-principles QCD description of heavy hadrons is still lacking due to the 
unknown nonperturbative QCD dynamics. 

In this study, the p-wave heavy meson is explored in a light-front quark model (LFQM) 
with both general and heavy quark frameworks. LFQM is a promising analytic method 
for solving the nonperturbative problems of hadron physics [34|], as well as offering much 
insights into the internal structures of bound states. The basic ingredient in LFQM is the 
relativistic hadron wave function which generalizes distribution amplitudes by including 
transverse momentum distributions; it contains all the information of a hadron from its 
constituents. The hadronic quantities are represented by the overlap of wave functions and 
can be derived in principle. The light-front wave function is manifestly a Lorentz- invariant, 
expressed in terms of internal momentum fraction variables which are independent of the 
total hadron momentum. Moreover, the fully relativistic treatment of quark spins and 
center-of-mass motion can be carried out using the so-called Melosh rotation [35[. This 
treatment has been successfully applied to calculate phenomenologically many important 



meson decay constants and hadronic form factors |36M40|| . Therefore, the main purpose of 
this study is the calculation of the leading twist LCDAs of p-wave heavy mesons within 
LFQM. 

The remainder of this paper is organized as follows. In Sec. II, the leading twist LCDAs of 
p-wave meson states are shown in cases of vector and tensor currents with general and heavy 
quark frames. In Sec. Ill, the formulism of LFQM is reviewed briefly; then, the leading 
twist LCDAs are extracted within the LFQM. In Sec. IV, numerical results of the decay 
constants and LCDAs are recorded. The £ moments of these LCDAs are also calculated and 
presented. Finally, conclusions are given in Sec. V. 

II. LEADING TWIST LCDAS OF p-WAVE MESONS 
A. General Framework 

Amplitudes of hard processes involving p-wave mesons can be described by the matrix 
elements of gauge-invariant nonlocal operators, which are sandwiched between the vacuum 
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and the meson states: 



(0\q(x)r[x,-x}q(-x)\H(P,e)), 



(2.1) 



where P is the meson momentum, e is the polarization vector or tensor (e does not exist in 
the case of the scalar meson), T is a generic notation for the Dirac matrix structure, and 
the path-ordered gauge factor is 



[x, y]=P exp 



ig, 



I dt(x - y)^{tx + (1 - t)y) 
Jo 



(2.2) 



This factor is equal to unity in the light-cone gauge which is equivalent to the fixed-point 
gauge, (x — y)^A^(x — y) = 0, as the quark- ant iquark pair is at the lightlike separation [41 ] . 
For simplicity, the gauge factor will not be shown below. 

The asymptotic expansion of exclusive amplitudes, in powers with large momentum trans- 
fer, is governed by the expanding amplitude Eq. (12. ip . shown in powers of deviation from 
the light-cone x 2 = 0. There are two lightlike vectors, p and z, which can be introduced by 



P =0, 



0. 



(2.3) 



such that p — > P in the limit M\ — > and z — > x for x 2 = 0. From this, it follows that J4] 

1 



x^ -P" 



M 2 H 



Px 



{Px) 



x 2 M 2 H 



x^ - P" 



X 



2Pz 



+ 0(x 4 ), 



2Pz' 



(2.4) 



where Px = P ■ x and Pz — pz — y (Px) 2 — x 2 Mfj. In addition, if it is assumed that the 
meson moves in a positive e% direction, then p + and z~ are the only nonzero components 
of p and z, respectively, in an infinite momentum frame. For the axial- vector meson, the 
polarization vector e M is decomposed into longitudinal and transverse projections as 



rf = « ( Jf - ,Mk' 



^ — f » _ f ^ 

1 e ± — e e || j 



pz \ 2pz 
respectively. For the tensor meson, the polarization tensor is 



(2.5) 



e^(m) = (11; m 1 W|ll; 2m) e»(m')e v (m"), 



(2.6) 
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(m is the magnetic quantum number) or 



4* = 4l4l: (2-7) 

4x = ^[4x^ + ^41], (2-8) 



'i^i+^ii + vl^' (2 - 9) 



and e /i .(= e^ u z u ) can also be decomposed into longitudinal and transverse projections as 

ejf* = — ^-^77^ , eT = e"* - ef. (2.10) 

LCDAs are defined in terms of the matrix element of a nonlocal operator in Eq. ( 12. ip . 
For scalar (S), axial vector (A), and tensor (T) mesons, the leading twist LCDAs can be 
defined as 



(0\q(z)Yq(-z)\S(P)} = f s / due^ z 



i 



o 



M 2 
2pz 



(2.11) 



(0\q(z)^ l5 q(-z)\A(P,e x=0 )} = if A M A f du e* pz {p»—<p A \\(u) + e» ±9A± (u) 

Jo 1 pz 



ez 

2 



—M 2 A g A3 (u)}, (2.12) 



(0\q(z)a^ l5 q(-z)\A(P,e x=±1 )) = [ du e^{(e^ - e^)4> A± (u) 



+ (p^-p» z »)^h Al{ (u) + (e^-elzn^h A3 (u)}, 

(2.13) 

(0\q(z)^q(-z)\T(P,e x=o )) = / T M| ['du ^{jT-^L^u) + € -^9t±(u) 

Jo 1 »zr 



" rM 2PF M|to(n) }' (2 - 14) 



(0\q(z)a^q(-z)\T(P,e x=±1 ))=tf^M T f du j^l^—^-fa^u 

Jo 1 pz 



(2.15) 



where w is the momentum fraction and £ = (1 —u)—u = 1 —2m. Here 05, (/)a,t\\ and 0a,t_l are 
the leading twist-2 LCDAs, and the others contain contributions from higher-twist operators. 
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The leading twist LCDAs are normalized as 

du0s,A(u) = l, (2-16) 

du^(j) T {u) = 1 (2.17) 

and can be parametrized as the so-called £ moments, 

(0 = t diCm- (2-18) 



1-1 

To disentangle the twist- 2 LCDAs from higher- twist operators in Eqs. ( 12. lip ~ (I2.15p . 
a twist-2 contribution of the relevant nonlocal operator q(z)Yq(— z) must be derived. In 
the case of Y — 7^(75), the leading twist-2 contribution contains contributions from the 
operators which are fully symmetric in Lorentz indices 42|, |43j : 

U-zmiMz)}2 = E rr?(o) i — -f Y + \ ,\ D»t ( 75 )g(o), (2.19) 

where D = ^ - ^ and 15 = — igB a (X a /2). The sum can be represented in terms of a 
nonlocal operator, 

[g(-^(75)^)]2 = Cdt^-q{-tz) fa 5 )q(tz). (2.20) 
Taking the matrix element between the vacuum and the p-wave meson state, we obtain 

(0\[q(-z)^q(z)] 2 \S(P)) = f s £ du<P s {u)\p^e^ z + (P" - p") £ dte* tpz W2I) 

(0|[g(-^)7 M 75g(^)] 2 |A( J P,e A= o)) = ^/iM, f l du4> A \\{u){v' i -e^ z 

Jo I 



+ ^ - p^^j £ dte* tpz X, (2.22) 



(0\lq(-z)Yq(z)h\T(P,e x=0 )) = f T M 2 T J q ^0 T|| ( M )|^__ e ^ 

+ 2 pV^5 / rfte^ p2 :2.23) 

ypz \PZ) J Jo 

For the derivations in Eqs. (I2.2ip ~ (I2.23p . we refer to Ref. 42], which dealt with the vector 
meson state. We can use Eq. ( 12.191) . and then expand the right-hand sides of Eqs. ( I2.2ip ~ 
(I2T23D . as follows: 

£ ^(0|g(0){^-^ + <Z ' ^ D* t) q mS{P)) 
Z^ n\ I n + 1 n + 1 J 
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/sEl du<p s {u)(&z) n U + (P^-jf) / dtt n , (2.24) 
n=0 n\Jo { Jo J 

E ^(0|#){^r7^ + w( ^ff" 1 gv}75g(0)|A(P, e )) 
^ n! I n + 1 n + 1 J 

tf A M A jr 1 ^ C du<j> A \ { {u){&z) n {^- + (e»-p>>—) f'dttA, (2.25) 
n! 7o [ ?2 \ pz J Jo J 

£l(0|gl0){^^ 7 - + ^^pDvW0)|T(P, e )) 
Z^ n! L n + 1 n + 1 J 

A*£S j£ ^<»)«*r{^ + 2 (£ - f t <«•}, (2.26) 



oo -n 

E 

n=0 

respectively. Picking n = in Eqs. f!2.24p and ( I2.25p . we obtain 



(0|g(0) 7 ^(0)|S(P)> = f s P» fdutsiu), (2.27) 

Jo 

(0|g(0)7 M 7 5 g(0)|A(P, eA=0 )) = i/^e" /' du<j> A \\{u). (2.28) 



Note that the tensor meson cannot be produced by the V — A current. We then pick n — 1 
in Eq. (I2.26P and obtain 



I(0|g(0)( 7 ^-P+ ^/>)g(0)|T(P, eA=0 )> = f T M^' £ du^ n (u). (2.29) 



From the normalization Eq. (I2.16D. w e have (0|g7 A1 75g| 3 y4i(P, e)) = ifz Al M3 Al e^ which is 



consistent with the results of Ref. 



M\- 



Next, we consider the case of V = (7^(75), where the leading twist-2 contribution contains 
contributions from the operators: 



V \ n+ \ P>V^]( 75 )g(0). (2.30) 



The sum can also be represented in terms of nonlocal operators: 



[q(-z)a^( l5 )q(z)h = dt 







^-q(-t 2 z)a^( l5 )q(t 2 z) + z a ^-q(-t 2 z)a^( l5 )q(t 2 ^ 



2.31) 



Taking the matrix element between the vacuum and the axial-vector and tensor meson state, 
we obtain: 

(0\{q(-z)a^ l5 q(z)] 2 \A(P,e x=±1 )) 
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el- 







/> 


4>A±(u) 







pz 



+ ( h A \\(u) - <f>A±.{u) 
(0\[q(-z)a^q(z)} 2 \T(P,e x = ±1 )) 



if^M T / duU T ±_{u) 



gVgi&z _|_ 



(2.32) 



pz 



+ [h T \\(u) - & T \{u 



pz 

T'^e^P* + 37*"" / dte^' 

V pz J Jo 



pz 



(2.33) 



where 



s ,pu 

rj-lp,V 



M 2 
2pz 



1 



2pz 



M 2 
pz 
M 2 



e" M 2 



u ,pv 



M 



pz 



^W-eV 4 ). 



(2.34) 



The derivations of Eqs. ( 12.321) and ( 12.33}) are shown in Appendix B of Ref . 45J . In contrast 
to Eqs. (I2.2ip ~ (I2.23p . the twist-2 LCDAs do not disentangle entirely from the higher 
twists in Eqs. (12. 32ft and (I2.33H . Taking the product with e± fJ ,z l , and e^,z u in Eqs. (12. 32ft 
and (I2.33p . respectively, we obtain 

1 



(0\[q(-z)a^e ± ^ 5 q(z)] 2 \A(P,e x=±1 )) = /J;jf du(f) A± (u)-(e • e ± Pz) \e* pz + jf dte** 



pz 



(0\[q(-z)^'e ± ,.q(z)] 2 \T(P,e x=±1 )) = i#M T / du<f> T± (u)-e^e ±tl . 



(2.35) 

e iipz + f 1 dte iit2pz 
Jo 

(2.36) 

Then, we use Eq. (I2.30P and expand the right-hand sides of Eqs. f !2.35|) and f !2.36|) as 

EA(°l^ ) (ra ^ ( "" D)n ^^75g(0)|A(P,6 A=±1 )) 



n=0 



711 



2n + l 



oo -n fl I [ rl 

= ZiE-T / du<p AL (u)-(e-e ± Pz)(Zpzr 1+ / 
£r n\ Jo 2 Jo 



dtt 



2n 



(2.37) 



00 I i 



n=0 



(n + l)(z-Dr ^, 

2n + 



oo -n r \ i f \ 

/^MrE-r/ ^0T±(«)-e M *^.(M n 1+ / ^ 2n 

n^O W ' "' 2 [ JO 



(2.38) 



S 



Picking n = and n = 1 in Eqs. ( 12 . 3 7|) and (12. 38 p . respectively, we obtain 

(0|g(0)^'e ±M75 g(0)|A(P,e A=±1 )) = fi t du<j> A± (u)(e ■ e ± Pz), (2.39) 

Jo 

(0\q(0)(z ■D)a^e ±lx .q(0)\T(P,e x=±1 )) = f^M T [' du^ TX (u)(e^e ±lx .Pz). (2.40) 

Jo 

B. Heavy Quark Framework 

In general, the theoretical description of meson properties relies on the bound state models 
with a relativistic normalization: 

(M(P', e')\M(P, e)> = 2P°(27r) 3 5 3 (P' - P)6 ee >. (2.41) 

At low energies, however, these models have little connection to the fundamental theory of 
QCD. The reliable predictions are often made based on symmetries. A well-known example 
is HQS [29], which arises since the Compton wavelength, 1/ttiq, of a heavy quark bound 
inside a hadron is much smaller than a typical hadronic distance (about 1 fm), and tuq is 
unimportant for the low-energy properties of the state. For a heavy-light meson system, it 
is more natural to use velocity v M instead of momentum variables. Then, it is appropriate 
to work with a mass-independent normalization of a heavy-light meson state: 

(M(v',e')\M(v,e)) = 2v°{2-k) 3 5 3 {Av' - Av)5gg>, (2.42) 

where A = M — itlq is the so-called residual center mass of a heavy-light meson. The relation 
between these two bound states is 



\M(P, e)) = VM\M(v, e)). (2.43) 



In addition, the heavy quark field can be expanded as [29] 
Q(x) = e - im Q v - x 



1 + i JPi 

iv ■ D + 2mq — is 



K(x), (2.44) 



where h*(x) is a field describing a heavy antiquark with velocity v. Then, the current qTQ 
can be represented as 

iTQ = '" r f 1 + S + ' ' ') h - (2 - 45) 
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Substituting Eqs. fl2^3|) and (12^51) into the definitions of LCDAs, Eqs. fl2TT|) ~ (l2~T5j) 
yield 



fOO 

F s I du e 1 ' 



(0\q(z)^h v (-z)\S(v)) 



(0\q(z)Y^h v (-z)\A(v,e x=0 )) = iF A / duo e 1 ' 

Jo 



2vz 

f z 



ez 



- z h 



2{vz) 



;G A3 (u) 



, (2.46) 



(2.47) 



(0\q(z)a^ l5 h v (-z)\A(v,e x=±1 )} — F A du j" vz (e>^ - e v ± v^ A± (u) 

Jo L 

ez 1 
- ^)— + (el^ - eT^)— i^aM 



(0|g(z) 7 ^„(-^)|T( V ,e A=0 )) = F T / e*' 

JO 



w vz 



2(vzf 

(0\q(z)a^h v (-z)\f(v,e x=±1 )) = iF^ I du e 



(2.48) 



(2.49) 



(iTv u - eW} 



vz 



+(// - v v z»)—,H n (uj) + (e^V - ?Z< 



[vz 



2{vz) 



H T3 (u)) 



(2.50) 



where Fm = vM/m, 3>i(w) = 4>i(u)/M, and u; was first introduced in Ref. [46j] as the 
product of the longitudinal momentum fraction u of the light (anti) quark and the mass of 
heavy meson M, namely, u = uM. Following a similar process, the leading twist LCDAs 
are obtained as 



(0^(0)7^(0)1^^)) = F s v» / du<S>s(tu 



(0|g(0) 7 M 75^(0)|i(t;,e A= o)) = iF A e» \ dcu<5> All (cu 



-{Q\q(0)(Yz-D+ ^)h v (0)\T(v,e x=0 )) = F T e^ J du$ n (u), 

(0|g(0)^*e ±M 75/i„(0)|^,e A=±1 )) = F^(e-e ± vz) / du<$> AX _{u 



(0\q(0)(z ■ D)a^'e ±IM .h v (0)\T(v,e x=±1 )) = F^e ±IM / du$ T± (uj). 

J 



o 

oo 



(2.51) 
(2.52) 
(2.53) 
(2.54) 
(2.55) 
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III. GENERAL FORMULISM IN LFQM 



A. General Framework 

A meson bound state, consisting of a quark qi and an antiquark q 2 with total momentum 
P and spin J, can be written as (see, for example, Ref. {37]) 



\M(P, L, J)) = J {d 3 h}{d 3 k 2 } 2(2tt) 3 5 3 (P - h - k 2 ) 



x 



X; n J s(h, k 2 , A 1; A 2 ) IqtxC^x, AO&Cfca, A 2 )), (3.1) 



Ai,A2 

where &i and &2 are the on-mass-shell light-front momenta, 



and 



~k = (k+,k ± ), fc ± = (fc 1 ,fc 2 ), fc~ = *\ \ (3.2) 



{d d k} : 



2(2tt) 3 ' 

|?(*i,Ai)g(A%, l A a )> =6l 1 (fc 1 )4 a (fc 2 )|0), (3.3) 
{b y (k'),b{(k)} = {d y (k'),d[(k)} = 2(2vr) 3 5\k' - k) 5 X 'X- 

In terms of the light-front relative momentum variables (u, Kj_) are defined by 

k+ = (1 - A;+ = «P+, 

ki± = (1 - u)P± + k±, k 2 ± = uP± - k±. (3.4) 
The momentum-space wave function ^ J L J S Z for a 2S+1 Lj meson can be expressed as 

^ J L f{k 1 ,k,X 1 ,X 2 ) = ^(LS;L z S z \LS-,JJ z )R s x ^ 2 (u,k ± ) <p LL ,(u, K± ), (3.5) 



where <^ll z (m, describes the momentum distribution of the constituent quarks in the 
bound state with the orbital angular momentum L, (LS; L Z S Z \LS; JJ Z ) as the corresponding 
Clebsch-Gordan coefficient and R\ t L constructing a state of definite spin (S, S z ) out of light- 
front helicity (Ai, X 2 ) eigenstates. Explicitly, 

R 1i\ 2 (. u > ^J-) = X ( a iI 7 ^m( 1 ~ M > Kx,mi)\s 1 )(X 2 \'R t M (u, -K ± ,m 2 )\s 2 )(- ~; sis 2 |~; SS$.6) 
where |sj) are the usual Pauli spinors, and T^m is the Melosh transformation operator 



3: 



{s\K M {u i ,K ± ,m i )\X} = — , (3.7) 

(rm + MiM ) 2 + k\ 
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with ui = 1 — Uj U2 = and n = (0, 0, 1) is a unit vector in the z-direction. In addition, 
Ml = ei + e 2 2 = ^ + -2 =t 

Ml « 2 

— r - ; — ; ei - k 2 e 2 + k z 

d = V m i + K ± + K *> — ; — = 1 - — ; — = w, 
v " ei + e 2 ei + e 2 

where k z is the relative momentum in the z direction and can be written as 

Km = ^ - ^t; 1 . (3.8) 

Mo is the invariant mass of qq and is generally different from the mass M of a meson which 
satisfies M 2 = P 2 . This is due to the fact that the meson, quark and antiquark cannot be 
simultaneously on-shell. We normalize the meson state as 

(M(P', J', J' Z )\M(P, J, J z )) = 2(2n) 3 P + 6 3 (P' - P)5 JIJ 5j /zJz , (3.9) 

in order that 

2(2?!-) ^ Lp '* L ' L '^ u ' k ^ll z (u, = h',L h' z ,L z - (3.10) 

Explicitly, we have 

VlL z = KL z Pp, (3.11) 



where kl z =±i = T( K ±x ±i/c_Ly) / '\/2, ^l z =o = K z are proportional to the spherical harmonics 
Yi£ z in momentum space, and is the distribution amplitude of the p-wave meson. In 
general, for any function F(\k\), <f p (u, n±) has the form of 



<P p (u,k ± ) = N^F(\k\), (3.12) 

where 

dKz = e ± e 2 . . 

du u(l - m)M 1 ' J 

is the Jacobian of transformation from (u, k±) to k and the normalization factor N is de- 
termined from Eq. (I3.10p . 

a p-wave meson state, it is more convenient to use the covariant form of 



T)SS Z 



In the case o 

nn 



36|,|39j,|47|: 

(IS; L Z S Z \1S; JJ Z ) k Lz R^ 2 (u, k ± 



SS Z / V K l K 2 



v/2 M (M + mi + m 2 ) 

xu^x, Ai)(f + Af )IWi P _v(fc 2 , A 2 ), (3.14) 
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where 

M = y/Mg - K - ™ 2 ) 2 - P = &i + fc 2 , 

A)w(fc, A') = ^^a,a', A) = ^ + m , 

«(*, A)«(Ar, A') = -^ A , A ,, £ A M fc ' A ) = ^~T+~~' (3-15) 

For the scalar, axial-vector, and tensor mesons, we have 

r.p, = £^7"(-if"), (3.16) 
where K = {k 2 — ki)/2 and 



£ A=±1 



— el(±l) • P ± , 0, e ± (±l) 



2 
.P 

e ± (±l) = T(l,±0/V2, 

6A =°- iWo V ^ J' ( } 

Note that the polarization tensor of a tensor meson satisfies the relations e^ u = e UfJi and 
e^P^ = = 0. Equations ( 13. 14ft and ( 13.161) can be further reduced by the applications of 
equations of motion on spinors: 



/ k k 

(lS;L z S z \lS;JJ z )k Lz Rf%(u,K ± ) = ^^u^, Ai)r^ +lp ^(A; 2) A 2 ), (3.18) 
where 



2 



3Po " 2v/3M : 
T' 1Pi = e • K75, 



IV = (« 2 - 26 • /f(rm - ma) j 7 5, 

/ 2^ \ 

r 3p 2 = ^ 7^ + U I ~ I ~ (-#")■ ( 3 - 19 ) 



M + mi + m 2 , 

Next, the matrix elements of Eqs. (1X27]) . (12T28D . (12T29|) . (T2~39l) . and (1X40|) will be calcu- 
lated within the LFQM, and the relevant leading twist LCDAs are extracted. For the scalar 
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meson state, we substitute Eqs. (13. ip . (13.51) . and (I3.18P into Eq. (I2.27P to obtain 



(0\frfqi\S(P)) = N c j{d 3 h} £ ^is z (ki,k 2 ,X 1 ,X 2 ){0\q 2 rq 1 \q 1 q 2 ) 

Ai,A2 



k + k + 
fsP* 1 I ducf>(u). 



r 



fa + rnA 



M 2 



fa + m 2 



kf J 2V3M 



(3.20) 



For the "good" component, \x = +, the leading twist LCDA <fts can be extracted as 

y/2W c f d 2 K ± [lM7li - (1 - w)m 2 ]M 



<t>s\u) 



(3.21) 



/s •/ 2(2tt) 3 ,^3n(l - u)M 

A similar process can be used for the axial-vector and tensor mesons which correspond 
to Eqs. (J22SD, (123511 . and (ESHD, fl23DJ, respectively, and the leading twist LCDAs are 
extracted as 



(p3 Al \\(u) 



f d 2 K± <p P (u, K±) 



M, 



<f)3 Al ±(u) 

(j>T\\{u) 



M J 2(2tt)3 ^ u (i_ u)Mo m L2M ( 

— 2« 2 (mi — m 2 )(umi + (1 — u)m 2 ) 
2^6 y d 2 k_l umi + (1 — u)m 2 



fi Al J 2(2tt) 3 ^(i _ u)M 

\/3 /" c/ 2 Kj_ ^(w, 



^-{M 2 - 4 mi m 2 ) + M 2 /^(l - 2u) 
o 

(3.22) 
(3.23) 



K p 



/4 ■/ 2(2vr)3 ^/^TT^Mo 

/" d 2 Kj_ K 2 , 



K i 



M (umi - (1 - tt)m 2 ) + 4(mi - m 2 )(3.24) 

M M 



/ii 7 2 (2tt) 3 0i(l - u)M Q 
v^6 /" (i 2 K^ fp{u, K±) 



<P P {U,K±), 



(3.25) 



h J 2(2tt) 3 ^(i _ U )M { 

+ (2k 2 z -k 2 ± )((1-2u) + 2 
2\/6 y (i 2 «:j_ <£>p(w, 



(ami — (1 — u)m 2 ) 
M + mi + m 2 



(3.26) 



/t 7 2(2tt) 3 _ u ) Mq 
2k 2 



K1 



K z (um 1 + (1 - w)m 2 ) - 77T-p(mi - m 2 ) 



+ 



M + m 1 + m 2 



(3.27) 
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B. Heavy Quark Framework 



If one takes mi = mg — > oo, that is, the heavy quark limit in the heavy-light meson, then 
two inequalities, uiq ~ M m 2 , n± and u — > 0, are obtained. The exact form of $ M can 
be derived by the redefinition of the meson bound state. Let us consider the bound states 
of heavy mesons in the heavy quark limit: 



\M(v;L,J)) = J{d 3 q}{d 3 k 2 }2(2n) 3 5 3 (Av-q-k 2 ) 

x ]T ^ J L f(u,K ± ,X 1 ,X 2 )bt(q,X 1 )d\k 2 ,X 2 )\0), 



(3.28) 



where q — ki — mqv is the residual momentum of the heavy quark. Operators bl(q, X±) 
create a heavy quark with 

{b v (q,X l ), bl,(q', X[)} = 2(2n)% v ,5 3 (q-q')5 XlK . (3.29) 

The variable u> and the relative transverse and longitudinal momenta, k± and k z , are ob- 
tained by 

u m\ + k\ 

u = e 2 + k z , n± = k 2± - ujv ± , k z = — — — . (3.30) 

The momentum-space wave function ^> J L J S Z can be expressed as 



1 



*^*(o;,« ± ,Ai,A 2 ) = -^(LS-LzS^LS^JJ^RlXiu,^) (3-31) 



where 

(LS; L Z S Z \LS, JJ z )k Lz R ss *(u, k ± , \ u X 2 ) = 
with 



V2^J (u + m 2 ) 2 + k] 



u(u,Ai)ru(fc2,A2) (3.32) 



f^Po = ~^(v ■ k 2 + m 2 ), 

fi Pl = 6 ■ /C 2 75, 

f 



3Pl - "71 

f 3 p 2 = -e^kz, 



v ■ k 2 + m 2 ) %-i-k 2 



75, 



(3.33) 



and 



Ml _ 
e A=0 — 



-e ± {±l)-v ± , 0, el(±l) 



-1 + v] 



(3.34) 
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u(v, Ai) is the spinor for the heavy quark, 

^u(v,\)u(v,\) = ^. (3.35) 

A V 

The normalization of the heavy meson bound states can then be given by 

(M(v', J', J' z )\M(v, J, J z )) = 2(2n) 3 v + 5 3 (Kv' - Av)5jj,8j z j,, (3.36) 

which not only leads to Eq. (12.43 j) . but also to the space part (pll z (u, k 2 ^) (called the light- 
front wave function) in Eq. (13.281) which has the following wave-function normalization 
condition: 



dud ,k± 



v?ll>,4)I 2 = i, ( 3 - 37 ) 



2(2tt) 3 

where (pih z = ^l z ^P p - In principle, the heavy quark dynamics are completely described by 
HQET, which is given by the 1/rriQ expansion of the heavy quark QCD Lagrangian: 

OO / 1 \ n 

£ = Q(iP- m Q )Q = £ -— C n . (3.38) 

n=0 V m Qj 

Therefore, \M(v; L, J)) and (p p (u, are determined by the leading Lagrangian Co = h v iv ■ 
Dh v . From the normalization conditions of Eqs. (13.101) and (I3.37p . we obtain the relation 
between wave functions <p p (u, k\) and <f p (u, k 2 ± ): 

ip p (u, = VM(fp(uj, k 2 ^). (3.39) 

In addition, in the heavy quark limit (mi — > oo), the heavy quark spin sq decouples from the 
other degrees of freedom so that sq and the total angular momentum of the light antiquark 
j are separately good quantum numbers. Hence, it is more convenient to use the L j j = P 2 3/2 , 



3/21/21/2 1 

P 1 , Pi and P basis. It is obvious that the first and the last of these states are 6 P 2 



and 3 P , respectively, while 



3; 



42 1 



i^ 3/! > = yfr^+^A). (3.40) 



In terms of P^ 2 and P^ 2 states, the relevant vertex functions read 



f p i/2 = -y=(v ■ k 2 + m 2 ) ^75, 

f p 3/ 2 = -—f=[(v k 2 + m 2 ) ^- 3e • /c 2 ]7 5 . (3.41) 
16 



Next, the matrix elements of Eqs. fl23T]) . fl23^ . fl233l) . fFCTj) and ([2755]) can be calculated, 
and the relevant leading twist LCDAs are extracted as 

F s ®s(u) = F 4 i/ a $ A i/2 || (w) = F^ /a $ i4 i/2_ L (w) 

= ^1 2(2% ~ 'J" 2 ^ + m2Y + K2±(pp{u ' ^ (3 - 42) 
F A w$ A w Ju) = f ^^^d^=L=[uj 2 -ml- 34 + 2k\ + 2ma«J, (3.43) 

^3/2 $.3/2 = / , f — [2^ - k ± \, (3.44) 



F T $ r ||(w) = v 7 ^ J 



d 2 K L <p p (u, k\) \ (u + m 2 ) 2 + K? ± 2 2 



(3-45) 



2 (2vr) 3 ^/(o; + m 2) 2 + ft:i . 

F^ T± ( W ) = V6 / f ^ > 2 - m 2 - 24 + m 2 «J, (3.46) 

7 2(2vr) 3 /( a; + m2 )2 + ft .2 



where we have denoted the P and P x states by Ay and A/ , respectively. In the 



literature, we read that 
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50j HQS requires 

f A v* = fs, / A 3/ 2 = 0. (3.47) 

These relations can be understood from the fact that (P^ 2 , Py^ 2 ) and (Py^ 2 , P^ 2 ) form two 
doublets in the heavy quark limit and that the tensor meson cannot be induced from the 
V — A current. From Eq. ( 13.421) . it is easy to find that 

F s = F A i/2 = Fj- 1/2 , $sM = ^v a „H = $ A i/ 2± (w). (3.48) 



Combining the former of Eq. (13.481) and the scaling F M = \/Mf M , we obtain: f s = f .1/2 = 

f\/2, which are consistent with the former of Eq. ( I3.47p . On the other hand, if wave 
A i 

function (p v has a similar form to Eq. ( 13 .121) . 



p (co,k ± ) = N'^F(\k\) } (3.49) 

where N' is the normalization constant and F(|k|) is a function of \k\, then, Eqs. (13.431) and 
(I3.44[) can be shown as F.3/2 = F\ /2 = or / 3/2 = f\ /2 = 0, which are also consistent with 

A l A x A l A 1 

the latter of Eq. ( I3.47p . The derivations are shown in Appendix A. In addition, the tensor 

meson can be created through the V — A currents with covariant derivatives [see Eqs. (12. 29ft 

and ( 12.401) ] . Thus, we can study its decay constant Ft and F^r here. From Appendix A, we 

find that Eqs. (ET35jl and (137361) lead to F T = F^ = and then ft = ft = f 3/2 = ft 3/2 = 0, 

1 A i 

which is consistent with the fact that (Py^ 2 , P^ 2 ) forms a doublet in the heavy quark limit. 
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IV. NUMERICAL RESULTS AND DISCUSSIONS 



In this section, the decay constants and LCDAs for the p-wave states of D, D s , B and 
B s systems are studied with the wave function <p p (u,K±). In principle, (p p is obtained by 
solving the light-front QCD bound-state equation Hlf\M) = M\M), which is the familiar 
Schrodinger equation in ordinary quantum mechanics, and Hlf is the light-front Hamil- 
tonian. However, except in some simple cases, achieving the full solution has remained 
a challenge. There are several popular phenomenological light-front momentum distribu- 
tion amplitudes which have been employed to describe various hadronic structures in the 
literature. A widely used one is the Gaussian type. If 

ISI2> 



F 9 (\k\) = exp 
then the corresponding wave functions are 



2/3 s 



(4.1) 



IY / 7T \ 3 / 4 



eie 2 



exp 



k\ + 



' uMp 
> 2 



2uM a 



2(3 2 



(4.2) 



y u(l - u)M 

Prior to numerical calculations, the parameters mi, 777.2, and /3, which appeared in the wave 
function, have to first be determined. Here, we use the parameters obtained in the ISGW2 



model 511 ] . an update of the ISGW quark model [52| for semileptonic meson decays. They 



are listed in Table I. Next, we use the parameters in Table I to evaluate the relevant decay 



TABLE I: The input parameters mi, mi, and (5 (in units of GeV) in the Gaussian- type wave 
function Eq. 



m u>d 


m s 


m c 


m b 


ficu 


Pes 


Pbu 


Pbs 


0.33 


0.55 


1.82 


5.20 


0.33 


0.38 


0.35 


0.41 



constants, Eqs. (13.2ip ~ (13.27p . for the p-wave heavy mesons. Note that the decay constant 

of the scalar and axial vector mesons are calculated by normalization of Eq. ( 12. 16ft . and 

that of the tensor mesons are calculated by normalization of Eq. (12. 17ft . The results, which 

compare with other theoretical evaluations, are listed in Table II. We find two things: the 

first is /^i/2 ^ /^3/2 ? an d the second is fa and /jr for the charm sector are larger than those 

for the bottom sector. Both are qualitatively consistent with the results in Section 3B, which 

were derived from the HQS. In addition, in this work, the ratio of the decay constants f.1/2 

A i 
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TABLE II: Mesonic decay constants (in units of MeV) obtained in this work and other theoretical 



evaluations. 



2S+1 Lj{L j j) 


3 p 

-t 


3 Pi 


l Pi 


3 P 2 


pl/2 


3/2 

1 


feu 


78 


-113 


45 


-50 


118 


-29 


feu [39J 


86 


-127 


45 




130 


-36 


feu [50] 


139 ± 30 








251 ±37 


77 ± 18 


feu [53] 


133 


-211 


72 








ftu 




-127 


80 


-48 


150 


-8 


fcs 


66 


-123 


38 


-65 


123 


-40 


fes [39J 


71 


-121 


38 




122 


-38 


fcs [50] 


110 ±18 








233 ± 31 


87 ± 19 


fcs [53J 


112 


-219 


62 








fis 




-107 


87 


-62 


138 


9 


fbu 


73 


-73 


42 


-17 


84 


-8 


fbu [39] 


112 


-123 


68 




140 


-15 


fbu [50] 


162 ± 24 








206 ± 29 


32 ± 10 


fbu [53] 


145 


-150 


76 








fbu 




-119 


52 


-17 


127 


-26 


fbs 


75 


-82 


43 


-25 


92 


-12 


fbs [5pJ 


146 ± 19 








196 ± 26 


36 ± 10 


fbs [53J 


140 


-157 


76 








fbs 




-122 


58 


-25 


133 


-23 



and fs in the D s system is 



/,y 2 //s = 1.9. 



(4.3) 



This is very close to the predictions of Refs. 5GJ and 54j: f A \/i/ fs — 2-12 and 2.26 ± 0.41, 
which are based on the mock-meson approach and the factorization hypothesis, respectively. 

For these heavy mesons, it is convenient to study the leading twist LCDAs using the 
Lj basis. Thus, the curves of </>(£) for the p-wave states in D, D s , B, and B s systems 
are evaluated and shown in Figs. 1 ~ 7. From these figures, we find the curves of 
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3 
2 



-1 

-2 

-1.0 -0.5 0.0 0.5 1.0 

FIG. 1: Leading twist- 2 LCDAs (f>s(0 of the heavy meson. The solid and dotted lines, and short 
and long dashes correspond to the D, D s , B, and B s systems, respectively. 

4 
3 
2 


-1 
-2 

-1.0 -0.5 0.0 0.5 1.0 

FIG. 2: Leading twist-2 LCDAs (ft .i/2,,(0 of the heavy meson. The solid and dotted lines, and 

A ! II 

short and long dashes correspond to the D, D s , B, and B s systems, respectively. 

A i/2||(£) and <fi A i/2 ± (£) are very similar to those of 0s(£), but are quite different from those 
of A 3/2||(£) and (j) A 3/2 ± (£). The scales of Figs. 3, 5, 6, and 7 are much larger than those of 
the others because the values of f^j 2 an d fx are relatively small. Finally, we parametrize 
the LCDAs in terms of the first six £ moments with Eq. (j2.18p . The results are shown 
in Tables III, IV, V, and VI. We find that the similarities between the longitudinal and 
transverse projections of the tensor meson are displayed not only in ratio fr/ fr — 1> but 
also in the approximations (C)t\\ — (C)t± f° r all heavy meson systems. 
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FIG. 3: Leading twist- 2 LCDAs (f).3/2,,(0 of the heavy meson. The solid and dotted lines, and 

A ! II 

short and long dashes correspond to the D, D s , B, and B s systems, respectively. 




FIG. 4: Leading twist-2 LCDAs (j> 1/2 (£) of the heavy meson. The solid and dotted lines, and 
short and long dashes correspond to the D, D s , B, and B s systems, respectively. 



V. CONCLUSIONS 



This study discussed the leading twist LCDAs of p-wave heavy mesons within the light- 
front approach. These LCDAs have been displayed in terms of light-front variables (u, u, k±) 
and the relevant decay constants in both general and heavy quark frameworks. In the 
heavy quark framework, we analytically found that the decay constants and LCDAs had 
the following relations: f s = f A i/2 = f^ 1/2 , $ s (w) = $ A i/^(w) = $ A i/2 ± (w) and f T = 
/jr = f A 3/ 2 = = 0, which are consistent with the requirements of HQS-(P 1 ^ 2 , P^ 2 ) 

and (P^ 2 , P^ 2 ) form two doublets. It was worth noting that we could study because 
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FIG. 6: Leading twist- 2 LCDAs (f>T\\{Q of the heavy meson. The solid and dotted lines, and short 
and long dashes correspond to the D, D s , B, and B s systems, respectively. 

the tensor meson was created through the V — A currents with covariant derivatives. In 
the general framework, we quoted the parameters mi and (3, which appear in the Gaussian- 
type wave functions, from the ISGW2 model, and numerically found that: i) an inequality 
/^i/2 3> /^3/2 is existent for all systems; ii) the decay constants fx and for the charm 
sector are larger than those for the bottom sector; iii) the ratio f.yz/fs = 1.9 for the D s 



system is close to the predictions in Ref. [50( and [54| : iv) the curves of ^i/ay ± (£) were very 
similar to those of 4>s(£), but were quite different from those of 4>j^/ 2 \\ an d 0t||,_l(6)> 
and v) the ratio Jt/Jt — 1 an d approximations (C)t\\ — (C)t± were satisfied for all heavy 
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20 
15 
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<frA& ^ 


-5 
-10 

-1.0 -0.5 0.0 0.5 1.0 

FIG. 7: Leading twist-2 LCDAs (j>T±(? of the heavy meson. The solid and dotted lines, and short 
and long dashes correspond to the D, D s , B, and B s systems, respectively. 





TABLE III: 


The first four 


£ moments of (f> 


M {0 for the 


p-wave states of D 


system. 




M 


S 


Pl /2 \\ 


pl/2 ± 


Pl /2 \\ 


p3/2 ± 


T|| 


T _L 


(? 

(?) 
(?) 
(?) 
(?) 
(?) 


0.14 
0.0029 
-0.050 
-0.061 
-0.062 
-0.057 


0.25 
0.091 
0.023 
-0.0031 
-0.015 
-0.019 


0.30 
0.13 
0.055 
0.023 
0.0061 
-0.0021 


1.7 
1.3 
1.0 

0.79 
0.62 
0.50 


-1.7 
-1.5 
-1.3 
-1.1 
-0.87 
-0.72 


1 

0.83 
0.68 
0.54 
0.43 
0.35 


1 

0.83 
0.68 
0.54 
0.44 
0.35 




TABLE IV: The first four £ moments of 4>m(Q for the p-wave states of D 


3 system. 




M 


S 


p l A\ 




p3/ 2|| 


p3/2 ± 


T\\ 


T _L 


(0 

(?) 
(?) 
(?) 
(?) 
(?) 


-0.13 
-0.14 
-0.14 
-0.11 
-0.092 
-0.072 


0.093 
0.011 
-0.029 
-0.032 
-0.030 
-0.026 


0.13 
0.033 
-0.013 
-0.020 
-0.022 
-0.020 


1.5 
0.96 
0.69 
0.48 
0.35 
0.26 


2.9 
1.9 
1.4 
1.0 
0.74 
0.55 


1 

0.66 
0.50 
0.35 
0.26 
0.20 


1 

0.66 
0.50 
0.36 
0.26 
0.20 



meson systems. It is easily realized that the results of i, ii, and iv qualitatively supported 
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TABLE V: The first four £ moments of 0m(£) for the p-wave states of B system. 



M 


s 


P} /2 \\ 
r i II 


pl/2 ± 


p3/ 2|| 
r l II 


p3/2 


T | 


T ± 


\s/ 


0.62 


0.63 


0.67 


2.4 


0.50 


1 


1 


(£ ) 


0.39 


0.41 


U.47 


2.9 


U.22 


1.4 


1 A 

1.4 


<£ 3 > 


0.25 


0.27 


0.33 


3.0 


0.059 


1.6 


1.6 




0.16 


0.18 


0.24 


2.9 


-0.034 


1.6 


1.6 


<e 5 > 


0.095 


0.12 


0.18 


2.8 


-0.088 


1.5 


1.5 


<e 6 > 


0.054 


0.077 


0.13 


2.6 


-0.12 


1.4 


1.4 



TABLE VI: The first six £ moments of 4>m(£,) for the p-wave states of B s system. 



M 




Pl /2 \\ 


pl/2 ± 


p3/2,| 

r \ II 


p3/2 L 


T|| 


T _L 


(0 


0.50 


0.54 


0.58 


2.3 


0.27 


1 


1 


<a 


0.25 


0.29 


0.35 


2.6 


-0.057 


1.3 


1.3 


<a 


0.11 


0.15 


0.22 


2.5 


-0.20 


1.3 


1.3 


(e) 


0.035 


0.077 


0.14 


2.3 


-0.26 


1.3 


1.3 


<a 


-0.0066 


0.032 


0.085 


2.0 


-0.28 


1.1 


1.1 


(e 6 ) 


-0.029 


0.0046 


0.053 


1.8 


-0.27 


1.0 


1.0 



the requirements of HQS. Due to the lack of relevant experimental data, the consistencies 
between our estimations and the predictions of HQS are important. 
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Appendix A: Some useful identities 



We consider an integration as 

dud 2 K^ <Pp(w, k 2 ± ) 



/ 9^3 7T , 9 = g(m 2 ,K z ,K ± ). (Al) 
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Substituting Eq. (13.491) with Eq. (1Alj) . we obtain 



„ T , f dud 2 K± \dK z F(\k,\) , 9x , „ x 

7 2(2tt)3 V du; / (u; + m2)2 + Ki yi ±J 1 J 



Taking the heavy quark limit for Eq. (I3.13p . we obtain y^f = y 77; an d Eq. (1A2j) can be 
rewritten as 

iV' / — — -J Vl y =g(m 2 ,K z ,K 2 ,). A3 

From Eq. (13.301) . the variables w, e 2 , and k z have the following relations: 

^2 "I - / . .v 

u = e 2 + k x , = e 2 - k z . (A4) 

Thus, Eq. f|A3[) can be rewritten as 

J 2(2tt) 3 / 2e2 (e 2 + m 2 ) 



Besides the function g, the only variable in Eq. ( 1A5j) is |/c| because e 2 = ym 2 + |k| 2 . 
Therefore, if we designate g to some specific function in Eq. flAll) . the integration can be 
made aware by symmetry. For example, if g = k z , 

dud 2 K ± <p P (u,Kj) 

2(2tt) 3 ^ (w + m2) 2 + K 2 Z ' 

because k 2 is an odd function. The second case is g = k 2 — k±/2, 

^ 2(2tt)3 ^( w + m2 )2 + K 2 ^ * 2 ; 1 ; 

because the contributions of k 2 ,, k 2 ^, and equal. The third case is g = 2K 2 (e 2 + K z ) = 

lo 2 — m 2 — 

G^ 2 /t ± y P (^,K±) / 2 _ 2_ 2 )= f dud 2 K ± ip p (u } K\) 2 

2(2tt) 3 ^( w + m2 )2 + K 2^ m2 K±j J 2(2vr) 3 ^ + ^2 + ^^' 1 j 

where Eqs. (IA6I) and (IA7[) are applied. We can employ Eqs. ( 1A6I) . (IA7I) and (IA8I) to prove 
that the integrations of u> for Eqs. (I3.43p . (I3.44p . (13 .45 p . and (I3.46P are all equal to zero. 
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